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Abstract 

Using the PCAC relation, we derive a compact formula for the pion decay con¬ 
stant F n in the nonlocal chiral quark model. For practical calculations this formula 
may be used both in the Minkowski and in the Euclidean space. For the pion mo¬ 
mentum —> 0 it reduces to the well known expression derived earlier by other 
authors. Using a generalized dipole Ansatz for the momentum dependence of the 
constituent quark mass in the Minkowski space, we express in terms of a single 
integral over the quark momentum fraction u. We interpret the integrand as a pion 
distribution amplitude 4>(u). We discuss its properties and compare with the 7rDA’s 
obtained in other models. 


1 Introduction 

Recent data from CLEO (1, and E791|2| experiments triggered a new wave of theoretical 
studies of the leading twist pion distribution amplitude pDA). On one side the data 
have been reanalyzed taking into account NLO perturbative QCD effects, as well as 
nonperturbative effects parameterized within the QCD light-cone sum rules On the 

other hand nonperturbative models u usi and lattice QCD ra-ra have been employed 
to calculate the 7 tDA from the relatively nonrestrictive physical assumptions. Here the 
dual nature of the pion, being the quark-antiquark bound state and the Goldstone boson 
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of the broken chiral symmetry at the same time, makes such calculations interesting by 
itself, even if the data is not yet decisive enough to distinguish between different models. 

Pion distribution amplitude is usually defined by means of the following matrix element 
(see e.g. [Ej): 


4>A U ) 



^L e -ir{2u-l){nP) 

71 


(0 \d(7ir)ifi'y 5 u (—nr) 17r + (P)) 


( 1 ) 


in the light cone kinematics where two quarks separated by the light cone distance z = 
2r along the direction n = (1,0, 0,-1) are moving along the light cone direction n = 
(1, 0, 0,1) parallel to the total momentum P. Here F n = 93 MeV. In this kinematical 
frame any four vector v can be decomposed as: 


V 

F = —n» + —v? + 

2 2 x 


( 2 ) 


with v + = n ■ v, v = n ■ v, and the scalar product of two four vectors reads: 


1 + - 1 - + - - 
V ■ W = —V W + —V W — V± ■ W±. 


(3) 


In Eq. © the path ordered exponential of the gluon field, required by the gauge invariance, 
has been omitted since we shall be working in the effective quark model where the gluon 
fields have been integrated out. 

In the local limit matrix element m reduces to 


<0 \A“(x)\ 7r fe (P)) = -iP^FJ ab e~ iPx 


(4) 


where 

H“(x) =7/5(a:)7 M 7 5 y^(^) (5) 

is the properly normalized axial vector current. 

In Refs. ITS] (p n (u) has been calculated in the effective chiral quark model in which 
quarks interact nonlocally with an external meson field 

f7^( s ) = e H5TV(i)/F, (g) 


and acquire a momentum dependent constituent mass 

M(k) = M k = MF[k) 2 . (7) 

M is a constituent quark mass of the order of 350 MeV and F(k ) is a momentum de¬ 
pendent function such that E(0) = 1 and F(k 2 —> oo) —> 0. Function F(k ) embodies 
nonperturbative effects due to the nontrivial structure of the QCD vacuum. Indeed, F(k ) 
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has been explicitly derived within the instanton model P2jJ. In Refs. |13| where the calcu¬ 
lations were performed in the Minkowski space (instanton model is inevitably formulated 
in the Euclidean metric) a convenient Ansatz for F(k) was used: 


F(k) = 


-A 2 


k 2 - A 2 + i 


ie 


n 


( 8 ) 


With this Ansatz (j) n (u) as well as higher twist 7rDA’s were calculated in Refs.J^j and j20| 
respectively. 

The problem is, however, that in the model with the nonlocal interaction (and momen¬ 
tum dependent quark mass M^) the axial current © does not exhibit PC AC m-m 
More drastically, a naive vector current 

(9) 

is not conserved. In order to restore these properties extra currents have to be added 
to A a and f22],[233- These new pieces modify both model expressions for F n and for 
(j) n (u). While the formula for F n is well known in terms of the Euclidean integral |23j.|25j: 


p 2= iv f d " k g Ml - 

“J (27T)- (kl + MlY 


( 10 ) 


(here M’ k = dMk/dk 2 ) the form of the wave function has been a subject of different studies 
with, however, contradictory results. For example the distribution amplitude obtained in 
Ref.|5 is very close to the asymptotic form 


CW=6«(i-«) (11) 

where u = k + /P + is the momentum fraction carried by the quark, whereas in ReE. h! if 71 
4>i t{u) = 1 . 

In the present work we derive the Minkowski space formula for F 2 for the modified 
axial current replacing the naive current in Eq.fjl}. Our formula, when continued to the 
Euclidean space, reduces to Eq. CUD- However, when evaluated in the Minkowski space 
by methods developed in Refs.|13|. it can be represented as an integral over du from 
an integrand which we interpret as 0„.( u ). This function does not resemble m and is 
compatible rather with the constant wave function of Refs. MU than with the result 
obtained in the same model e I2], however, with the naive current e. 

There are several comments which are due at this point. First of all it is not clear 
how the modified current can be generalized to the bilocal operator entering formula o. 
That is why it was argued in Refs. ME that rather than considering matrix elements 
of the form m or ®D, one should calculate the whole physical process in the effective 
model, impose Bjorken limit to make contact with the expressions known from QCD and 
extract the distribution amplitude. One has to note however, that the effective models 
are not valid at large momenta which are needed to impose Bjorken limit. Moreover it is 
not clear whether the distribution amplitudes defined that way are universal. Secondly, 
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arguments may be given that it is not necessary to insist that the bilocals defining the 
distribution amplitudes must reduce to the proper currents in the local limit 1 * . Indeed, as 
we shall show below the naive bilocal m reproduces the Pagels-Stokar formula (2(1 for 


Ft: 


Ft = 4 Ah 


d 4 k E Ml - \k\M k M' k 
(2ir)‘ ( k% + Mlf 


( 12 ) 


which was obtained from the Ward-Takahshi identities. 


2 Currents in the nonlocal models 

Let us consider the model defined by an action ra.ra: 

S= j(dk)-f(k)(f;-m)il/(k)-M j(dkdl)ii(l)F{l)lPHl - k)F(k)ii(k). 

Here, following |^j ( dk ) = d 4 h/(27r) 4 etc., and ( dx ) = d A x. Equations of motion for the 
quark fields read 

— M [ ( dl)F{k)U l5 (k - l)F(l)i/j(l) + mx/j(k), 


= M J (dl)^(l)F{l)W 5 (l - k)F(k) + rmp(k). 
To get the equation of motion for the U 75 field let us expand © 

U l5 (l — k) — (2tt) 4 5 (4) (/ -k) + ^7 5 t c tt c (/ -k) + ... 

-*7T 

and the equation of motion gives a constraint 

[ ( dl ) ${l + k)F{l + k)-f 5 r a F (1)^(1) = 0. 


(13) 


(14) 

J (15) 

It is easy to verify that the naive vector current © is not conserved |23‘, |HjJ • I 11 order 
to restore current conservation, the following two currents have to be added to V[t (in 
momentum space) 

v;(P) = v;(p) + p;(p) + l%p) ( 16 ) 

with left and right currents dehned as 

Z 

L a (P)=iM [(dxdydz ) [ d Sfl e iPs $(x)F{x - z)T a U 15 (z)F(z - y)ip(y), 


R a (P ) = iM I ( dxdydz ) / ds^e iPs ip(x)F(x - z)U l5 (z)T a F(z - y)^(y) (17) 


1 By local limit we understand the limt in which the fields in Eq.(l) are taken in the same point x. 

There are still corrections due to the momentum dependent constituent mass and nonlocal interactions. 


4 




where T a = r“/2. Accordingly the modified axial current reads: 


A%{P) = A^P) + BRP) - L%{P) (18) 

with T a = , y 5 T a /2. The integral ds M should be understood as an integral over the path 
connecting points z and y or x. This prescription makes the L“ and currents path 

dependent [233 (strictly speaking the transverse part is not fixed). 

The divergence of the modified vector current is, however, path independent and takes 
the following form: 


P^VVP) = m / ( dkdl)ip(k)F(k ) 


r 


—, U l5 (k — l + P) 


mm- 


(19) 


This is immediately zero for the baryon current (r“ = 1). For the isospin current we can 
expand £/ 75 m 


T 




= —^T a e abc 7T c (k -l + P) + ... 

Fir 


and (USD vanishes due to the constraint (©. For the axial current we get 

p^A a jp) = -m [ mmrnrn + p ) 


-m i (dkdi)m)F{k)m—m{k~i+p)) F (i)m- ( 2i ) 


( 20 ) 


By expanding U 75 m we arrive at 


p^A a (p) = -m / mm) mm + p) 


m j (dk)m) F i k ) mp(k + p)m + p) 

M 


F v 


(dkdl)m)F{k)F(l)m ir a (k -l + P) + ... 


( 22 ) 


which is the proper PCAC formula (note that the second term vanishes due to (fTHIP . 

In order to calculate F n we can either use Eq.(|2J) with A a —> A“ |2Zj or use the PCAC 
relation 

(23) 


0 


id^A a Ax) 7r 6 (P) ) = -iP z FJ ab e 


ra6 —iPx 


which is what we are going to do in this work. Notice, that we have to calculate the 
matrix element in Eq. P31> off-shell, extract the leading power in P 2 and take the limit 
P 2 -»• 0. 
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Figure 1: Diagrams contributing to the matrix element of Eq. (21). Black squares denote r a,b 75. 

3 Decay constant and the distribution amplitude 

3.1 Matrix elements 


There are three contributions to the matrix element of Eq. (12 lj) depicted in FigQ one 
from the first term of expansion m and two (which by the anticommutation rule reduce 
to one term, see Eq. El) from the term in m involving one pion field. Adding all of 
them we get 


id^A^z) 


AP) 


<-qfe -iPz 

F 

- L 7T 


X 


(24) 



M u k MML - k M ± MiMizZ , Ml 

. ‘ t - P (B-Ml)((k-Pf-Ml_ p ) + B-Ml 


Symmetrizing the last term with respect to the change of variables k —> k — P , adding all 
terms and comparing with Ea. (I23H we arrive at 


Ft 




- MM - MizFAL 

[v - MD({k - py-- Mi_ P y 


(25) 


By expanding Eq. El in powers of P 2 we recover the Minkowski version of Eq. m- 
Indeed, by changing the variables: k —> k + P/2 we get 


F 2 , = -t4 N C P j(ik)\ 


P w fc+P/2^ 




2 ((k + £)* - M| +p/2 )((fc - £y - Mi_ Pn ) 


(26) 


In fact an expression identical to Eq. El appears in the axial and pseudoscalar corcllators 
derived within the instanton model of the QCD vacuum [25j . 

Noting that 

M k ±p/2 = M k ± ( kP)M' k + ~^M' k + ... (27) 

where ' denotes d/dk 2 we have 


1 




P 2 M 2 - A(Pk) 2 M k M' k + Ak 2 (kP) 2 Ad k 


/ 2 


M 2 k+P/2 


m 


Pi 

2 > 


M k ~p/ 2 ) 


(28) 
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Since under the integral k fl k l/ —> \g, u ,k 2 (plus a term proportional to P^P U which we may 
safely neglect) equation (0H1) transforms into 


Fl = -i4N c / (dk) 


Ml - k 2 M k M' k + k 4 Ml 2 


= —iAN c / (dk) 


((* + f ) 2 - ^ + P/2)((fc - f ) 2 - M fc-P/ 2 ) 

M 2 - k 2 M k M' k + k 4 M' k 2 


j (* 2 - M 2 ) 2 

On the other hand, matrix element of the naive current © , gives m 

Mk-pkfj, + M fc (P M - kfj) 


FlP„. = —iAN c 



kM k - P 


(P - Ml) ((* - P? - Ml_ p ) 


(29) 


(30) 


which by the same steps which led from Eq. (©3) to flu reduces equation (EHD to the 
Minkowski version of the Pagels-Stokar formula m- 


3.2 Calculation of the loop integral 

In order to calculate the loop integral in Eq. 425D with M k given by Eqs. (EE! > we shall 
introduce the light-cone parameterization of the momenta © with 

d 4 k = —^-dudk~ d 2 kp (31) 

where k + = uP + . The method of evaluating dk~ integral, taking the full care of the 
momentum mass dependence, has been given in 1©|. To evaluate dk~ integral we have to 
find the poles in the complex k~ plane. It is important to note that the poles come only 
from the momentum dependence in the denominators of Eas. (l25l30D . This means that 
the position of the poles is given by the zeros of denominator, that is by the solutions of 
the equation 

ta - M2 G»-£ +fc r +fc= °- (32 > 

This equation is equivalent to 

4n+l 

G(z) = z 4n+1 + z 4n - r 2 = \J(z- Zi), (33) 

i— 1 

with ^ = k 2 / A 2 — 1 + ie and r 2 = M 2 / A 2 . For r 2 ^ 0 (or finite A) equation i|33l) has 4n +1 
nondegenerate solutions which we denote Zi. Equation E2D should be understood as an 
equation for k~ = k~(zi). In general case 4 n of s can be complex and the care must be 
taken about the integration contour in the complex k~ plane. Because of the imaginary 
part of the Zi s, the poles in the complex k~ plane can drift across Re/c - axis. In this 
case the contour has to be modified in such a way that the poles are not allowed to cross 
it. This follows from the analyticity of the integrals in the A parameter and ensures the 
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vanishing of 7 tDA’s in the kinematically forbidden regions. The results are expressed as 
sums over zf s which have to be found numerically. 

In order to avoid spurious divergences coming from k~ in the numerator of Eq. © 
we shall make use of the Lorentz invariance, writing 


Fl = -i4 N c —I^\ 

Since I^ v vanishes for P M —> 0 (see (ESI) and ©) we have that 


with 

Then 


I,, = A(P 2 )P,P U + -B( y P 2 )P 2 g lll 


A(P 2 ) -> A, B(P 2 ) -> B for P 2 -> 0. 


F 2 = -i4N c (A + B ) 


Hence we have to calculate 2 integrals: 


A = j^n ,J n /l I fliy , B = V 


The result reads 

A = 


4n+l 


167T 2 


M 2 duY fifk(zl n u + z 2n u) 2 x In (1 + z t u + z k u ), 

0 i,k =1 

1 


P = 


i 2M 2 A 2 
16tt 2 P 2 


x In ( (1 + ztu + z k u) - —tm 


Here u = 1 — u and 


4n+l ^ 4n+l ^ 

/.=n—=n— 

1 1 * - ** l Zk ~ z * 

k^i k^i 


for which the following identities hold 

4n+l 


£/■ 


z m = 


i =1 


0 m < An 


1 rn — An 


(34) 

(35) 

(36) 

(37) 

(38) 

(39) 


4,1+1 / p2 \ 

du Y fifk{4 n - 4 n ) 2 ( (! + z iU + Z k u) - —uu J 
i,k =1 ' ' 


(40) 


(41) 


As seen from Eq. flUD the first term in B is singular as P 2 —> 0 in apparent contradiction 
with the finiteness of P 2 . However, the function 


N r M 


2 4n+l 


0inf(w) =- A— Y fifk{4 n - z 2n ) 2 (l + ZiU + z k u ) ln(l + ZiU + z k u) (42) 


i,k =1 








vanishes when integrated over du. Hence the finite formula for F 2 reads 


n m 2 r 

F l = —^5- du Yl K z l nu + z i * u ) 2 - 2 ( z l n - z 2n ) 2 uu\ In (1 + ZiU + z k u). (43) 

Q l,k=l 

This allows us to define the distribution amplitude 
N M 2 4 " + ! 

M U ) = ~ 4 \p2 hf k + ^ 2n “) 2 ~ 2 ( Z k n ~ ^") 2 H 111 i 1 + + Z kU) • (44) 

n i,k =1 

Let us recall that the distribution amplitude defined by means of the naive axial current 
(J3) reads JJh 

AT ][/[2 

M u ) = hf k ( z ?zl n u + zf n zy) In (1 + Zi u + z k u ). (45) 

7r i,k 

3.3 Numerical results 
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Figure 2: a) dots: cutoff parameter A for different n (Eq.(8)) and for M = 350 MeV, dahed line: fit 
described in the text; b) function </>inf(w) for n = 3/2, 3 and 5. 

Condition dm, or equivalently © provide a relation between parameter A, con¬ 
stituent mass M and power n from Eq.(JBJ). Throughout this paper we shall use M = 350 
MeV. The value of parameter A = A(n) obtained from Ea. (I43D . or from Eci. (1101) after 
continuation of the cutoff formula (13 to the Euclidean metric, is depicted in FigJ2]a. It 
is interesting to note, that our formula © for F 2 , unlike equation ©, does allow for 
half integer n’s. An approximate relation, depicted by a dashed line in FigEla holds 

A [MeV] = 432.82 + 444.61 n - 28.02 n 2 . 
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The local current © contributes, through Eq. G2D, approximately 70% to the total nor¬ 
malization. 

Having fixed A for given n, we can calculate the distribution amplitude as defined 
by Eq. (HU). However, before doing this we have to check whether the formally divergent 
part, given as an integral over u from the function </>i n f(n), vanishes. We have checked 
numerically that this is indeed the case. Function </> in f(-u) is plotted in FigOlb for n = 3/2, 
3 and 5. 

Next, in Fig© we plot the distribution amplitude <j) n (u) for n = 3/2 and n — 5 
(solid lines) together with the contributions from integrals A and B 137(1 . We see that 
the contribution from A is relatively flat and does not vanish at the end points. The 
contribution from B vanishes at the end points and is even negative in their vicinity. 
There is not much difference between the two cases n — 3/2 and n — 5, although one may 
say that the smaller n the flatter <p n . 




Figure 3: Function <p n for n = 3/2 and n = 5 with two contributions A and B given by Eqs.(39). 

In Fig Ola we plot for comparison function ^(u) (OH), o-(u) 03) corresponding to 
the naive axial current © for M = 350 MeV and n — 3, together with the asymptotic 
distribution amplitude (OH). One should note that while model distributions are defined as 
some low normalization scale Q 2 = /r 2 , </>“ s (w) corresponds to the limit Q 2 —> oo. Indeed, 
the leading twist distribution amplitude can be expanded in terms of the Gegenbauer 
polynomials 


Q 2 ) = 6n(l 


u ) 


OO 

1+ XI a n (Q 2 )Cl /2 (2u - 1) 

n= 2,4... 


(46) 


where a n (Q 2 ) —> 0 in the large Q 2 limit [28] . It is important no notice that a n {Q 2 ) tend 
to zero monotonically, so that they cannot change the sign. 
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Figure 4: a): Functions (Eq.(44), solid), (j>„ (Eq.(45), dashed) for n = 3 and asymptotic distribution 
amplitude (Eq.(ll), dash-dotted); b): Evolustion of <j) v from the initial scale Q = 425 MeV (dotted line) 
to the CLEO point Q = 2.4 GeV (solid) and for Q = 10 12 GeV (dashed) for n = 1. 


As soon as we switch on the QCD evolution, distribution amplitude 4> n (u, Q 2 ) changes 
the shape and it goes to zero at the end points. This evolution is plotted in Fig0b 
for n — 1, assuming 2 light flavors, A qcd = 175 MeV, and initial scale p = 425 MeV. 
This initial scale has been adjusted in such a way, that the second Gegenbauer coefficient 
a 2 , when evolved to the CLEO point Q = 2.4 GeV, gives a 2 (2.4) = 0.15 as indicated 
by the analysis of RefdSJ- For the normalization scale as large as 10 12 GeV the evolved 
distribution is slowly approaching the asymptotic one. 

4 Summary and discussion 

In the present paper we have derived a compact formula fU) for the pion decay constant 
in the nonlocal chiral quark model. In order to define F v we have used the full nonlocal 
axial current m and the PCAC relation (H. Equation K25D, when expanded in the pion 
momentum P M , reduces to the well known j23|. |25j formula (1291) . The advantage of Eq. P5|) 
consists in the fact that it can be evaluated in the Minkowski space with a suitable Ansatz 
for the momentum dependence of the constituent mass (JHJ). By integrating (EH) over dk 
and d 2 k_ j_ we are left with a du integral over the function (j) n {u) which we interpret as a 
pion distribution amplitude (HU). 

As mentioned at the end of Sect |2 it is not clear how to extend the local current m 
to the bilocal operator like the one entering formula ©. Therefore our definition of the 
pion distribution amplitude may be not correct. However, it is worth to note that the 
shape of our distribution amplitude resembles a constant DA obtained by a consistent use 
of the Ward-Takahashi identities 0-P, rather than the DA calculated in the instanton 
model of the QCD vacuum in Ref. |5J, although in both approaches full nonlocal currents 
have been used. 


11 









Unfortunately, the 7 tDA derived here and in Refs. 0-EH is probably phenomenolog¬ 
ically unacceptable. That is because the detailed analysis of the CLEO data indicates 
that the coefficient a 4 (2.4 GeV) is negative |Hj,|[4J and possibly as large as a 2 (2.4 GeV) 
(4j. In our case, however, a 4 is always positive. The same concerns the constant 7rDA. In 
this respect 7rDA derived in by the same methods in Refs. ra using the bilocal operator 
m with no extra pieces corresponding to the nonlocal currents m fits the data much 
better. That is because, similarly to the results of Refs. [Sj , it exhibits a shallow minimum 
around u = 1/2 which generates negative a 4 . 

As already mentioned above, there is a problem how to define the distribution am¬ 
plitudes in the effective models of QCD. This is due to the fact that the QCD currents 
and the model currents are not the same. One way would be to perform factorization 
and large Q 2 expansion in QCD and then parameterize the nonperturbative matrix ele¬ 
ments by a set of unknown distribution amplitudes. To calculate these matrix elements 
an effective model, like the one discussed here, is used. Considering operators as obtained 
from QCD leads to the violation of PCAC and, in the worse case, to the violation of the 
gauge invariance at the level of the effective model. Another method consists in perform¬ 
ing factorization and large Q 2 expansion directly in the effective model. This is possible, 
since the degrees of freedom of the effective models discussed here are, at least as the 
quantum numbers are concerned, identical to the degrees of freedom of QCD (except for 
gluons, which are not present in the former case). This means, however, that the low 
energy model has to be applied to the processes with large momenum transfer. Since the 
currents of the effetive models are not the same as in QCD, extra pieces contributing to 
the DA’s, as compared to the previous method, are present. Although in this work we 
have not calculated the physical process and have not implemented the Bjorken limit, our 
approach is in our opinion equivalent, since we have considered the matrix element © 
of the full current m ■ Our results indicate that these two methods lead to completely 
different DA’s . The first method gives the 7 tDA resembling the asymptotic distribution, 
whereas the second approach generates the DA which is compatible with a constant. 


Acknowledgments 

It is a pleasure to dedicate this work to Jan Kwiecinski. 

We would like to thank A. Rostworowski for comments and for reading the notes. M.P. 
would like to thank W. Broniowski and E. Ruiz-Arriola for comments and discussion. 
Special thanks are due to K. Goeke and all members of Inst, of Theor. Phys. II at 
Ruhr-University where part of this work was completed. M.P. acknowledges support of 
the Polish State Committee for Scientific Research under grant 2 P03B 043 24. 


References 

[1] J. Gronberg (CLEO Collaboration), Phys. Rev. D57 (1998) 33 hep-ex/9707031 


12 



[2] E.M. Aitala et al. (Fermilab E791 Collaboration), Phys. Rev. Lett. 86 (2001) 4768 
hep-ex/00100431; D. Ashery, talk at Workshop on Physics with Electron Polarized 
Ion Collider, hep-ex/99100024. 

[3] A. Schmedding and O. Yakovlev, Phys. Rev. D62 (2000) 116002 jhep-ph/9905392 . 

[4] A.P. Bakulev, S.V. Mikhailov and N.G. Stefanis, Phys. Rev. D67 (2003) 074012 
|!hep-ph/0212250]; hep-ph/0303039 

[5] A.P. Bakulev, S.V. Mikhailov and N.G. Stefanis, Phys. Lett. B508 (2001) 279 
hep-ph/0103119|; in: Proceedings of the 36th Rencontres de Moriond, QCD and 
Hadronic Interactions , hep-ph/0104290, 

[6] T. Heinzl, Ligh-Cone Quantization: Foundations and Applications , Lect. Notes 
Phys. 572 (2001) 55 |hep-th/0008096 ; Nucl. Phys. Proc. Suppl. 90 (2000) 83 
jhep-ph /0008314]. 

[7] T. Shigetani, K. Suzuki, H. Toki, Phys. Lett. B308 (1993) 383 |hep-ph/9402286|; 
Nucl. Phys. A579 (1994) 413 hep-ph/9402277]. 

[8] A.E. Dorokhov and L. Tomio, Phys.Rev. D62 (2000) 014016; I.V. Anikin, A.E. 
Dorokhov and L.Tomio, Phys. Lett. B475 (2000) 361; Phys. Part. Nucl. 31 
(2000) 509, (Fiz. Elem. Chast. Atom.Yadra 31 (2000) 1023), A.E. Dorokhov, talk 
at 37th Rencontres de Moriond, QCD and High Energy Hadronic Interactions, 
hep-ph/0206088 

[9] R.M. Davidson and E. Ruiz Arriola, Acta Phys. Polon. B33 (2002) 1791 
hep-ph/0110291]; Phys.Lett. B348 (1995) 163. 

[10] E. Ruiz Arriola, Acta Phys. Polon. B33 (2002) 4443 hep-ph/0210007|. 

[11] E. Ruiz Arriola and W. Broniowski, Phys.Rev. D66 (2002) 094016 |hep-ph/0207266|; 
Phys.Rev. D67 (2003) 074021 hep-ph/0301202|. 

[12] V.Yu. Petrov and P.V. Pobylitsa, hep-ph/9712203; V.Yu. Petrov, M.V. Polyakov, R„ 
Ruskov, Ch. Weiss and K. Goeke, Phys. Rev. D59 (1999) 114018 |hep-ph/9807229|. 

[13] M. Praszalowicz and A. Rostworowski, Phys. Rev. D64 (2001) 074003 
hep-ph/0105188|; talks at 8th Adriatic Meeting, Particle Physics in the New Millen¬ 
nium, hep-ph/0202226 and at 37th Rencontres de Moriond, QCD and High Energy 
Hadronic Interactions, hep-ph/0205177, 

[14] D. Daniel, R. Gupta and D.G. Richards, Phys. Rev. D43 (1991) 3715. 

[15] L. Del Debbio et al, Nucl. Phys. Proc. Suppl. 83 (2000) 235 |hep-lat/990914Tj|. 


13 









[16] M. Burkardt and H. El-Khozondar, Phys. Rev. D60 (1999) 054504 |hep-pli/98054951; 
M. Burkardt and S. Seal, Phys. Rev. D65 (2002) 034501 hep-ph/0102245 ]; 
hep-ph/0101338 

[17] S. Dailey, Phys. Rev. D64 (2001) 036006 hep-ph/0101318|; 

[18] P. Ball, JHEP 9901 (1999) 010 jhep-ph/9812375 . 

[19] D.I. Diakonov and V.Yu. Petrov, Nucl. Phys. B245 (1984) 259; B272 (1986) 457. 

[20] M. Praszalowicz and A. Rostworowski, Phys.Rev. D66 (2002) 054002 
Ihep-ph/01111961; Acta Phys. Polon. B34 (2003) 2699 [hep-ph/0302269 . 

[21] B. Holdom, J. Terning and K. Verbeek, Phys. Lett. B232 (1989) 351, B. Holdom, 
Phys. Rev. D45 (1992) 2534. 

[22] G. Ripka and R. Ball, in: proceedings of Conference on Many-Body Physics, Coimbra 
1993, hep-ph/9312260 

[23] R.D. Bowler and M. Birse, Nucl. Phys. A582 (1995) 655 hep-ph/9407336|; R.S. 
Plant and M. Birse, Nucl. Phys. A628 (1998) 607 |hep-ph/9705372|. 

[24] W. Broniowski, in: Miniworkshop on Hadrons as Solitons, Bled 1999, 
hep-ph/9909438 

[25] D.I. Dyakonov and V.Yu. Petrov, Sov. Phys. JETP 62 (1985) 431 (Zhur. Eksp. Teor. 
Fiz. 89 (1985) 751). 

[26] H. Pagels and S. Stokar, Phys. Rev. D20 (1979) 2947. 

[27] A. Bzdak, in preparation. 

[28] A.V. Efremov and A.V. Radyushkin, Theor. Mat. Phys. 42 (1980) 97; Phys. Lett. 
B94 (1980) 245; S. Brodsky and G.P. Lepage, Phys. Lett. B87 (1979) 594; Phys. 
Rev. D22 (1980) 2157. 


14 





